function y = source(x)
    y = sin(pi*x);
end